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A  simple  model  is  employed  to  illustrate  the  influence  of  the  dynamic  effective  potential  resulting  from  torsional  motion  on 
tunneling  via  bending  modes  in  trans-cis  isomerization  in  polyatomic  molecules.  Wc  show  that  the  effect  is  significant  and  can 
determine  the  dominant  tunneling  mode  in  systems  such  as  HSiOH. 


1.  Introduction 

Most  treatmenis  of  barrier  tunneling  make  use  of 
simplifying  assumptions,  such  as  one-dimensional 
paths.  As  chemists. we  are  often  interested  in  pro¬ 
cesses  that  occur  in  polyatomic  molecules,  thus  the 
tunneling  processes  of  interest  are  often  multi-di¬ 
mensional.  Makri  and  Miller  [1  ]  have  presented  a 
method  for  treating  tunneling  in  multi-dimensional 
systems  semiclassically.  However,  it  is  often  neces¬ 
sary  and  appropriate  to  calculate  tunneling  correc¬ 
tions  by  assuming  that  it  occurs  along  some  special 
"tunneling  direction",  for  example,  along  a  normal 
mode  direction. 

An  example  where  this  approach  may  be  reason¬ 
able  is  trans-cis  isomerization.  It  is  often  a  good  ap¬ 
proximation  to  treat  the  torsional  motion  as  sepa¬ 
rable.  However,  as  ab  initio  calculations  have  shown, 
trans-cis  isomerization  can  occur  via  an  angle  bend¬ 
ing  mode  in  some  molecules.  The  purpose  of  this 
Letter  is  to  point  out  that  there  are  dynamic  effects 
which  arise  from  the  torsional  motion  that  must  be 
taken  into  account  for  tunneling  along  a  bending 
mode  direction.  Wc  use  a  very  simple  model  to  il¬ 
lustrate  the  qualitative  effects  of  the  torsional  me- 

1  Permanent  address:  Institute  of  Theoretical  Chemistry.  Jilin 
University.  Changchun,  1 .10023,  Peoples  Republic  of  China. 

;  Current  address:  Department  of  Chemistry.  University  of 
Oklahoma.  Norman,  OK,  USA. 


tion  on  the  bending  mode  tunneling.  We  consider  the 
simplest  case,  that  is.  a  four-atom  molecule,  namely. 
HSiOH. 

The  cis  and  trans  contormcrs  of  HSiOH.  and  us 
isotopic  (D  and  lsO)  analogs,  have  received  a  lot  of 
attention  lately.  Ismail  et  al.  [2]  studied  the  inter¬ 
action  of  silicon  atoms  and  w  ater  molecules  in  a  solid 
argon  matrix.  Infrared  spectra  of  the  products  (in¬ 
cluding  cis  and  trans  conformers  of  HSiOH. 
HSilsOH.  HSiOD.  and  HSi'sOD)  were  measured 
and  assigned.  Withnall  and  Andrews  [  3 )  studied  the 
matrix  reactions  of  silane  and  oxygen  atoms  and 
identified  HSiOH.  HStlsOH.  DSiOD.  and  DSilsOD 
by  infrared  absorption.  Kudo  and  Nagase  [4.5]  have 
reported  ab  initio  calculations  for  H;SiO  and  related 
molecules,  including  trans-  and  cis-HSiOH.  They 
found  that  there  are  two  possibilities  for  the  trans  to 
cis  isomerization  transition  state:  either  along  an  in¬ 
plane  bending  mode  or  along  the  torsional  direction. 
Their  results  showed  that  the  latter,  which  has  a  bar¬ 
rier  of  1 1  kcal/mol.  is  1.1  keal/mol  lower  than  the 
former.  Sakai  and  Jordan  [6]  have  reported  SCF  re¬ 
sults  for  the  siructures.  stabilities,  and  normal  mode 
frequencies  for  cis-  and  trans-HSiOH.  ( They  also  re¬ 
ported  results  for  similar  molecules.  They  found  that 
the  lowest  energy  path  for  trans -»cis  isomerization 
in  HAIOH  involves  the  in-planc  motion  of  the  H- 
atom  bonded  to  the  O-atom  rather  than  the  torsional 
motion  as  in  the  case  of  HBOH.  HCOH.  and 
HSiOH.)  Tachibana  et  al.  [7.8]  have  investigated 
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tunneling  in  the  cis-trans  isomerization  in  HSiOH 
be  using  an  ab  initio  IRC  (intrinsic  reaction  coor¬ 
dinate)  method.  They  proposed  that  the  isomeri¬ 
zation  path  may  be  regarded  as  one-dimensional 
along  the  direction  of  the  SiOH  bending  mode.  That 
is.  the  path  maintains  Cs  symmetry  and  the  v  ibra¬ 
tional  modes  orthogonal  to  it  do  not  contribute  to 
the  dynamics  along  the  path. 

We  should  also  point  out  some  classical  dynamics 
studies  that  are  relevant  to  this  problem.  Spears  and 
Hutchinson  [9. 10  j  used  classical  mechanics  to  study 
the  dynamics  of  vibrationally  excited  trans-ditmide. 
Csing  a  very  simple  model  and  classical  trajectory 
calculations,  they  showed  that  the  isomerization  dy¬ 
namics  of  diimide  are  dominated  by  the  "fictitious'' 
forces  arising  from  the  presence  of  torsional  motion 
and  the  source  of  the  final  "kick”  of  energy  into  the 
isomerization  coordinate  |  torsion )  arU^i  from  a 
Coriolis  coupling  of  the  torsion  to  the  NH  stretches. 
They  pointed  out  that  the  isomerization  via  a  planar 
bending  mechanism  is  prevented  since  such  motion 
entails  crossing  an  infinite  effective  potential  due  to 
the  centrifugal  force.  That  is.  they  concluded  that  the 
bending  pathway  is  inactive. 

We  focus  here  on  the  dy  namic  effects  on  tunneling 
paths  between  the  trans  and  cis  conformers  of  mol¬ 
ecules  of  the  type  HABC.  taking  HSiOH  as  an  ex¬ 
ample.  We  use  a  simple  model  to  illustrate  the  effect 
ot  the  dynamic  effective  potential,  which  arises  from 
the  torsional  motion,  on  the  tunneling  probability.  It 
is  pointed  out  that  this  "dynamic"  correction  to  the 
"static"  potential  in  one-dimensional  treatments  of 
tunneling  for  cis-trans  isomerizations  may  be  im¬ 
portant.  We  also  note  that  due  to  the  existence  of  this 
dynamic  effect,  the  realistic  tunneling  path  may  be 
a  curved  route  that  is  a  combination  of  the  torsional 
and  bending  modes,  and  is  possibly  coupled  with  the 
stretching  modes  as  well. 

2.  The  dynamic  effective  potential  barrier 

The  Born-Oppenheimer  approximation  allows  the 
Schrodinger  equation  for  a  molecule  to  be  separated 
into  two  equations. 


ir,  +  r,,+  i  ;.n  +  rnnm.(*.r) 

=  E(R)<PAlt.r)  .  (1) 

{Tn  +  l-(R)]'PjR)  =  eVjR)  .  (2) 

where  y/c(R.  r )  and  ip„(R)  are  the  electronic  and  nu¬ 
clear  wavefunetions.  respectively.  Rand  rare  the  po¬ 
sition  vectors  of  the  nuclei  and  electrons,  respec¬ 
tively:  Tc.  7„.  I  «..  1  cn.  and  I  „„  are  the  operators  for 
the  kinetic  energies  for  the  nuclei  and  electrons,  elec¬ 
tron  repulsions,  electron-nuclear  attractions,  and 
nuclear  repulsions,  respectively  :  c  is  the  total  energy 
and  /:(/?)  is  the  "potential  energy  surface"  (which 
we  refer  to  as  the  static  potential). 

Now.  let  us  consider  the  trans  -cis  isomerization 
in  HABC-type  molecules  taking  HSiOH  as  an  ex¬ 
ample.  The  ab  initio  calculations  [4,7]  show  that  the 
trans-  and  cis-HSiOH  conformers  exist  in  local  min¬ 
ima  with  planar  structures  on  the  ground-state  po¬ 
tential  energy  surface.  The  two  minima  are  sepa¬ 
rated  by  barriers  of  about  10  keal/mol  along  both 
the  out-of-plane  torsional  coordinate  and  the  in-plane 
bending  mode  direction.  Kudo  and  Nagase  [4]  and 
Sakai  and  Jordan  [6]  have  suggested  that  the  tor¬ 
sional  mode  is  the  dominant  pathway  for  the  iso¬ 
merization.  However.  Tachibana  et  al.  [7.8]  as¬ 
sumed  that  the  isomerization  involves  pure  SiOH 
bending  motion  and  computed  tunneling  probabil¬ 
ities  for  passage  through  the  barrier  along  that  co¬ 
ordinate.  The  dominani  pathway  for  the  isomeri¬ 
zation  is  determined  simply  by  the  relative  heights 
of  the  two  barriers  on  the  static  potential  energy  sur¬ 
face.  In  the  following,  we  iust  show  that  the  effective 
potential,  which  arises  from  dynamic  effects,  should 
be  considered  for  this  process. 

The  coordinate  system  for  our  model  is  illustrated 
in  fig.  I.  The  SiO  bond  is  along  the  c  axis  with  the 
O-aiom  at  the  origin.  Points  A.  B.  C.  D.  and  E  are 
all  located  on  a  sphere  of  radius  R(l.  the  equilibrium 
OH  bond  length.  The  equilibrium  cis  and  trans  ge¬ 
ometries  both  lie  in  the  )'OZ  plane.  The  OH  bond 
points  to  the  upper  right  for  the  trans  geometry  and 
the  upper  left  for  the  cis.  As  the  molecular  transla¬ 
tion  and  rotation  are  sepaiated  out.  the  static  poten- 
tia1  can  be  written  generally  as  I  '( r,.  r:.  0,.  I)2,  tp,  —  tp2 ). 
where  rt ,  fl,.  and  ip ,  describe  the  first  H-atom.  and  r2, 
02,  and  <p2  the  second.  The  torsional  potential  de¬ 
pends  only  on  the  difference  <pt  -<p2. 
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Fig.  I.  Iliuslration  of  the  coordinates  for  the  model.  All  atoms  are 
in  the  YOy  plane  at  the  cis  and  trans  equilibrium  geometries. 
The  points  A.  B,  C,  D.  and  E  represent  points  along  the  isomeri¬ 
zation  pathways,  where  A  and  C  are  the  trans  and  cis  conformers. 
respectively.  Path  A  -FI  ■(  is  the  pure  in-plane  bending  mode 
isomerization  Passage  through  point  E  involves  a  combination 
of  SiOH  bending  and  torsional  motion. 

The  torsional  motion  involves  rotation  of  the  two 
hydrogen  atoms  about  the  r  axis.  The  r  component 
of  the  total  internal  angular  momentum  l.  is 

T- =  -  .  (3) 

where  the  subscripts  1  and  2  indicate  the  two  hy¬ 
drogens.  The  potential  commutes  with  /... 

[  !  (/•,./•,.  0,.  02,  (/>,  -  ip2 )./..]  =0  .  (4) 

i  e.  L,  is  a  constant  of  the  motion  in  the  ground  state 
for  which  the  eigenvalue  of  L:  is  zero  since  the  total 
rotation  has  been  separated  out.  However,  the 
operator 

l-:  =  Ll+L;.  +  L:  (5) 

is  not  a  constant  of  the  motion,  that  is, 

[  H',.  r2.  0,.02.  (a,  -40,),  /.-]  *0  .  (6) 

From  symmetry  considerations  of  group  theory,  F( r,. 
r2.  0,.  th,  is  invariant  under  rotation  about 

the  r  axis  since  becomes  40,+ <p  and  <p:  becomes 
'o,  +  tp  as  the  system  rotates  by  an  angle  ip  about  the 
a  axis.  Thus.  V  should  commute  with  the  operator  /., 
since  the  infinitesimal  rotation  operator  about  the  r 
axis  is  I  —  i tL:.  where  t  is  the  infinitesimal  rotation 
angle.  Or  directly,  given  a  function  V(rt.  r2.  0,.  02, 


V>i-  <p2).  it  is  easy  to  verify  that  (VL:  —  L: 1  ).  since 
L.  =  —  lfid/dip,  —  \fid/d<p2  and  we  have 

/..(  '=().  From  this  it  is  clear  that  Lx  and  and  thus 
L -,  do  not  commute  with  the  potential  I  .  Thus,  the 
L 2  operator  has  a  distribution  of  values  rather  than 
a  definite  eigenvalue  for  the  internal  nuclear  motion 
of  HABC  in  the  ground  state. 

Now.  let  us  return  to  our  model  as  illustrated  in 
fig.  1.  Since  the  torsional  motion  (  r  in  fig.  I  )  is  an 
internal  rotation  about  the  SiO  bond  in  HStOH.  we 
can  regard  the  HSi  bond  fixed  in  the  FOZ  plane  and 
thus  consider  only  the  relative  rotation  of  the  H-atom 
bonded  to  O-atom.  When  this  H-atom  moves  from 
point  A  to  C  in  fig.  1 ,  the  geometry  changes  from  trans 
to  cis.  The  H-atom  moving  along  a  certain  path  s 
around  the  sphere  experiences  a  dynamic  effective 
potential  Fdyn  due  to  torsional  motion, 

ld>n  =  TJ/2p(»sinf/)- .  (7) 

where  n  is  the  mass  of  the  H-atom  and  the  other 
quantities  are  as  defined  in  fig.  1.  This  comes  from 
the  kinetic  energy  operator  7n  of  eq.  (2).  Thus,  the 
bending  motion  of  the  H-atom  is  coupled  to  the  tor¬ 
sional  motion  such  that  it  cannot  be  treated  as  purely 
one-dimensional  on  the  static  potential  energy 
surface. 

In  this  simple  model,  motion  along  the  reaction 
coordinate  s  and  torsional  direction  \  is  considered 
in  the  kinetic  energy  operator  T„.  Along  a  given  tun¬ 
neling  path  7n  =  t(p;+/^)  +  7',v  +  F(.v).  where 
j/>;  and  J/>;  are  kinetic  energy  operators  along  s  and 
v,  respectively.  7„  is  the  coupling  of  the  two  mo¬ 
tions.  and  F(.v)  is  the  potential  energy  along  the  v. 
Here,  we  only  want  to  explore  the  dy  namic  effect  due 
to  the  torsion  so  we  ignore  the  coupling  term  7's>.  The 
5 pi  term  is  represented  by  Z.’/Z.nfrsinf!)2.  where 
rsin  0  is  the  distance  from  the  H-atom  to  the  r  axis 
and  r  is  the  OH  bond  length.  In  more  accurate  treat¬ 
ments.  such  as  those  by  Marcus  [11]  and  Miller  et 
al.  [12],  the  reaction  path  Hamiltonian  can  be  writ¬ 
ten  (considering  only  the  coordinates  v  and  v)  as 


ni _ 

j  1  +Mv).v]- 


+  +  1 »(')  +  I  i  (s.  y)  . 


where  M  \ )  is  the  curv  ature  of  the  path  s.  v  is  the  dis¬ 
tance  deviated  normal  to  x.  F0(.v)  is  the  potential 
along  the  path,  and  I  ,  (  v.  v)  is  the  potential  change 
as  the  motion  deviates  from  the  path  v.  At  v=0.  i.e. 
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along  the  reaction  coordinate  v.  (', (v.  ,v)=0  and 
A ( v  ).y = 0.  which  leads  to  eq.  (7). 

It  should  be  noted  that  the  case  we  treat,  the  effect 
of  torsional  motion  on  the  tunneling  path,  is  differ¬ 
ent  to  those  usually  considered  [11.12].  The  usual 
treatment  is  for  vibrational  modes  that  are  normal 
to  the  path  v  with  a  harmonic  or  modified-harmonic 
model,  and  where  the  vibrational  energy  only  about 
the  path  s  is  considered. 

Let  us  now  make  a  rough  estimate  of  the  effect  of 
the  dynamic  effective  potential  on  the  tunneling  path 
in  the  trans~>cis  isomerization  of  HStOH.  We  refer 
again  to  fig.  1.  which  illustrates  isomerization  along 
the  pathways  with  torsional  motion.  When  the  hy¬ 
drogen  atom  is  at  point  A  r=0  and  0=65.5*  .  at  point 
B  r  is  undefined  and  0=0: .  at  point  C  r=  180  and 
tf=65.5  .  at  point  D  r  is  between  0  and  90  and 
0=  65.5  .  and  at  point  E  r=90  and  l)<  65.5  .  Pure 
planar  isomerization  is  along  the  path  ABC.  The 
more  general  path  involves  a  torsion  plus  bend  mo¬ 
tion  as  illustrated  by  the  points  ADEC  in  fig.  1.  where 
it  is  v  isualized  as  being  broken  down  into  simple  seg¬ 
ments  of  the  two  motions.  We  can  roughly  estimate 
this  path  by  assuming  it  follows  the  smallest  polar 
angle  0=  (65.5  /'90  )  r  as  the  atom  moves  from  A 
to  C.  Quantum  mechanically,  the  internal  rotation 
operator  I.2  has  an  average  deviation  (uncertainty* 
on  the  order  of  ft2.  Thus,  the  dependence  of  the  larg¬ 
est  value  of  the  dynamic  effective  potential  along  the 
path  on  r  can  be  expressed  approximately  as 

1  ;„P<  n  -  fi-^2/(/?(y  T)J 

=  0.05 sin  2(0.73r)  keal/mol .  (8) 

The  cis  and  trans  conformers  are  at  about  the  same 
energy  and  the  barriers  for  the  two  isomerization 
channels  on  the  static  potential  energy  surface  are 
roughly  at  the  same  level  (about  10  keal/mol  above 
the  conformer  wells).  We  approximate  the  barrier 
height  of  the  total  effective  potential  as 

r„(r)  =  [  10.0  +  0.05  sin  '2(0."’3r)  1  keal/mol  .  (9 ) 

in  order  to  simply  show  the  dynamic  effect.  The  de¬ 
pendence  of  (  ,,( r)  on  r  is  plotted  in  fig.  2.  It  de¬ 
creases  rapidly  as  the  torsional  angle  increases. 


Fig.  2.  The  barrier  height  of  the  approximate  effective  potential 
as  a  function  of  the  torsional  angle  r.  eq.  ( 9 ). 

3.  Tunneling  probability 

To  estimate  the  effect  of  the  dy  namic  effectiv  e  po¬ 
tential  on  the  tunneling,  we  represent  the  potential 
with  the  Eckart  barrier  [13]. 

l  (.s)=  I0(r)sech:(v)  .  <  1 0 ) 

where  v  is  the  tunneling  coordinate.  Fig.  3  shows  the 
shapes  of  the  potential  barrier  for  a  few  r.  where  the 
lowest  one  is  the  static  potential. 

The  one-dimensional  Schrodinger  equation  along 
the  tunneling  path. 

d2ii/(v)  2m  , 

-JV  +  t-,-  [/:  -r„  sech2(.v)  ]t</(  v)  =0  .  (  11  ) 

d.v  ■  fi  ■ 


s 


Fig  .3,  The  Eckart  barriers  used  to  represent  the  static  potential 
and  the  dynamic  effective  potentials  for  torsional  angles  r=  5  , 
10  .  15  .and  20  . 
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can  be  solved  analytically  [13].  The  permeability  is 
f,un=sinh:(  izk) 

X  [sinh:(  nk )  +  cosh-v'  1  jt(  &»il'0/fr  —  1 )  ]  ~  1  . 

(12) 

where 

k=x/2mE/f) . 

The  dependence  of  tunneling  probability  on 
the  torsional  angle  r  for  £=8  kcal/mol  is  shown  in 
fig.  4.  The  probability  of  tunneling  is  zero  at  r=0 
since  the  total  effective  potential  is  infinite.  As  r  in¬ 
creases  1  dcp  decreases  rapidly  (see  fig.  2'  and  the 
tunneling  probability  becomes  significant  at  about 
t=20\ 

The  torsional  energy  levels  are.  of  course,  quan¬ 
tized  with  relatively  large  spacings.  We  have  not  con¬ 
sidered  the  quantized  levels  here,  although  one  could 
easily  do  so.  We  merely  wish  to  illustrate  the  effect 
and  thus  have  considered  the  trans-cis  motion  on  a 
potential  energy  surface  where  the  torsional  mode 
might  be  considered  to  beany  state,  generally,  not  an 
eigenstate  of  the  torsion. 

We  hav  e  ignored  the  effects  of  the  variation  in  the 
width  of  the  barrier,  which  is  certainly  significant. 
We  have  assumed  that  all  paths  have  the  same  width 
since  we  are  interested  in  merely  illustrating  the 
qualitative  dynamic  effect  of  the  torsional  motion  on 
the  tunneling  paths.  Work  is  in  progress  on  a  more 
accurate  and  elaborate  treatment  of  this  problem. 


Fig  4  The  tunneling  probability  as  a  function  of  the  torsional 
angle. 


4.  Conclusions 

We  have  used  a  simple  model  to  illustrate  the  ef¬ 
fect  of  the  dynamic  effective  potential  that  results 
from  torsional  motion  on  tunneling  paths  along  the 
trans-*cis  isomerization  coordinate.  This  effect  may 
be  important  in  many  systems  and  must  be  explicitly 
considered  in  one-dimensional  treatments  of  tun¬ 
neling.  Because  of  this  effect,  the  torsional  or  a  curv  ed 
torsional-plus-bend  pathway,  possibly  coupled  to 
stretching,  is  likely  to  be  preferred  ov  er  the  pure  in¬ 
plane  bending  pathway. 
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